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Abstract 

We present a review of several results concerning the construction 
of the Cramer-von Mises and Kolmogorov-Smirnov type goodness- 
of-fit tests for continuous time processes. As the models we take 
a stochastic differential equation with small noise, ergodic diffusion 
process, Poisson process and self-exciting point processes. For every 
model we propose the tests which provide the asymptotic size a and 
discuss the behaviour of the power function under local alternatives. 
The results of numerical simulations of the tests are presented. 
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1 Introduction 



The goodness-of-fit tests play an important role in the classical mathematical 
statistics. Particularly, the tests of Cramer-von Mises, Kolmogorov-Smirnov 
and Chi-Squared are well studied and allow to verify the correspondence of 
the mathematical models to the observed data (see, for example, Durbin 
(1973) or Greenwood and Nikulin (1996)). The similar problem, of course, 
exists for the continuous time stochastic processes. The diffusion and Pois- 
son processes are widely used as mathematical models of many evolution 
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processes in Biology, Medicine, Physics, Financial Mathematics and in many 
others fields. For example, some theory can propose a difi^usion process 



dXt = S, {Xt) dt + a dWt, X^, 0<t<T 

as an appropriate model for description of the real data {Xt, < t < T} and 
we can try to construct an algorithm to verify if this model corresponds well 
to these data. The model here is totally defined by the trend coefficient 
5"* (•), which is supposed (if the theory is true) to be known. We do not 
discuss here the problem of verification if the process {Wt, < t < T} is 
Wiener. This problem is much more complicated and we suppose that the 
noise is white Gaussian. Therefore we have a basic hypothesis defined by the 
trend coefficient {■) and we have to test this hypothesis against any other 
alternative. Any other means that the observations come from stochastic 
differential equation 

dXt = S (Xt) dt + a dWt, Xq, 0<t<T, 

where S* (•) 7^ 5* (•). We propose some tests which are in some sense similar 
to the Cramer-von Mises and Kolmogorov-Smirnov tests. The advantage of 
classical tests is that they are distribution-free, i.e., the distribution of the 
underlying statistics do not depend on the basic model and this property 
allows to choose the universal thresholds, which can be used for all models. 

For example, if we observe n independent identically distributed random 
variables {Xi, . . . = with distribution function F (x) and the basic 
hypothesis is simple : F {x) = F^ {x) , then the Cramer-von Mises and 
Kolmogorov-Smirnov D„ statistics are 

-12 

F„ (x) - F^ (x) dF* (x) , = sup F„ (x) - F^ (x) 



respectively. Here 



1 " 



^ ^ {Xj<x} 



is the empirical distribution function. Let us denote by {Wq (s) , < s < 1} 
a Brownian bridge, i.e., a continuous Gaussian process with 

EWo (s) = 0, EWo (s) Wo{t)=tAs- st. 
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Then the hmit behaviour of these statistics can be described with the help 
of this process as follows 



sup |VFo(s) 

0<s<l 



Hence the corresponding Cramer-von Mises and Kolmogorov-Smirnov tests 
with constants Ca,da defined by the equations 



a, 



sup \Wo (s)| > dc 

Q<s<l 



a 



are of asymptotic size a. It is easy to see that these tests are distribution- 
free (the limit distributions do not depend of the function (■)) and are 
consistent against any fixed alternative (see, for example, Durbin (1973)). 

It is interesting to study these tests for nondegenerate set of alternatives, 
i.e., for alternatives with limit power function less than 1. It can be real- 
ized on the close nonparametric alternatives of the special form making this 
problem asymptotically equivalent to the signal in Gaussian noise problem. 
Let us put 

F {x) = (x) + ^ / h (F, (y)) dF, {y) , 
where the function /?,(•) describes the alternatives. We suppose that 

I /i(s)ds = 0, I h{sf As <oo. 
Jo Jo 

Then we have the following convergence (under fixed alternative, given by 
the function /i (•)): 



n 



h (v) dv + Wo (s) 



ds, 



sup 

0<s<l 



h (v) dv + Wo (s) 



We see that this problem is asymptotically equivalent to the following signal 
in Gaussian noise problem: 



dYs ^ K {s) ds + dWo {s) , < s < 1. 



(1) 
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Indeed, if we use the statistics 



r-l 



D — sup 



0<s<l 



then under hypothesis h{-) =0 and alternative /i (■) 7^ the distributions of 
these statistics coincide with the hmit distributions of and \fnDn under 
hypothesis and alternative respectively. 

Our goal is to see how such kind of tests can be constructed in the case 
of continuous time models of observation and particularly in the cases of 
some diffusion and point processes. We consider the diffusion processes with 
small noise, ergodic diffusion processes and Poisson process with Poisson and 
self-exciting alternatives. For the first two classes we just show how Cramer- 
von Mises and Kolmogorov-Smirnov - type tests can be realized using some 
known results and for the last models we discuss this problem in detail. 



2 Diffusion process with small noise 

Suppose that the observed process is the solution of the stochastic differential 
equation 



where Wt, < t < T is a Wiener process (see, for example, Liptser and 
Shiryayev (2001)). We assume that the function S{x) is two times contin- 
uously differentiablc with bounded derivatives. These are not the minimal 
conditions for the results presented below, but this assumption simplifies the 
exposition. We are interested in the statistical inference for this model in 
the asymptotics of small noise : £ — > 0. The statistical estimation theory 
(parametric and nonparametric) was developed in Kutoyants (1994). 

Recall that the stochastic process = {^^^^,0 <t<T} converges uni- 
formly in t G [0,T] to the deterministic function {xt,0 <t < T}, which is a 
solution of the ordinary differential equation 





(2) 



dxi 
lit 



t 



S{xt) 



0<t<T. 



(3) 
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Suppose that the function (x) > for x > and consider the following 
problem of hypotheses testing 



. 



S (x) = (x) , xo < X < x^ 
S (x) 7^ S** (x) , Xo < X < x^ 



where we denoted by x^ the solution of the equation ([3D under hypothesis 



Xo+ «) dv, 



0<t<T. 



Hence, we have a simple hypothesis against the composite alternative. 

The Cramer-von Mises {W^) and Kolmogorov-Smirnov [D^) type statis- 
tics for this model of observations can be 



w: 



dt 



S^xtY 



dt 



S4x;y 



^1/2 



sup 

0<t<T 



dt, 



It can be shown that these two statistics converge (as e 
functionals 



0) to the following 



w: 



W(sYds, e-^D 



sup |VF , 

0<s<l 



where ,0 < s < 1} is a Wiener process (see Kutoyants 1994). Hence 

the corresponding tests 

with the constants c^, da defined by the equations 



W(sY ds > Cn 



a. 



sup \W (s)| > da 

0<s<l 



a 



(4) 



are of asymptotic size a. Note that the choice of the thresholds Cq, and da 
does not depend on the hypothesis (distribution-free). This situation is quite 
close to the classical case mentioned above. 

It is easy to see that if S (x) ^ (x), then supo<t<T \xt — x\\ > and 



Wt 



oo, e — > oo. Hence these tests are consistent against any fixed 
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alternative. It is possible to study the power function of this test for local 
(contiguous) alternatives of the following form 

dXt = S,{Xt) dt + e^^dt + edWt, 0<t<T. 

We describe the alternatives with the help of the (unknown) function 
h (■). The case h{-) = corresponds to the hypothesis J^. One special class 
of such nonparametric alternatives for this model was studied in lacus and 
Kutoyants (2001). 

Let us introduce the composite (nonparametric) alternative 

Mi : h{-)enp, 

where 

Hp^i^hi-) : ^ \{xf //(dx) >p|. 

To choose alternative we have to precise the "natural for this problem" dis- 
tance described by the measure (■) and the rate of p = p^. We show that 
the choice 

[x) 

provides for the test statistic the following limit 

/■I r /"S "I 2 

^ / K (v) dv + W (s) ds, 
Jo Uo 

where we denoted 



ds 



K (s) = ul^'^h , = / 2 

Jo (x*) 

We see that this problem is asymptotically equivalent to the signal in white 
Gaussian noise problem: 

dYs = K (s) ds + dW (s) , 0<s<l, (5) 

with the Wiener process W {■). It is easy to see that even for fixed p > 
without further restrictions on the smoothness of the function (•) the 
uniformly good testing is impossible. For example, if we put 

hn {x) — cS^ {x)^ cos [n {x — Xq)] 
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then for the power function of the test we have 

inf </3(^s,/in) 
h{-)eHp 



a. 



The details can be found in Kutoyants (2006). The construction of the uni- 
formly consistent tests requires a different approach (see Ingster and Suslina 
(2003)). 

Note as well that if the diffusion process is 

dXt = S {Xt) dt + ea (Xj) dW*, Xq^xq, 0<t<T, 
then we can put 



T 



a [X 



dt 



\£S^{x*f 



dt 



and have the same results as above (see Kutoyants (2006)). 



3 Ergodic diffusion processes 

Suppose that the observed process is one dimensional diffusion process 

dXt = S (Xt) dt + dWt, Xo, 0<t<T, (6) 

where the trend coefficient S (x) satisfies the conditions of the existence and 
uniqueness of the solution of this equation and this solution has ergodic 
properties, i.e., there exists an invariant probability distribution Fs (x), and 
for any integrable w.r.t. this distribution function g (x) the law of large 
numbers holds 

fJ^gW dt^ P 9{x) dFs{x). 

These conditions can be found, for example, in Kutoyants (2004). 

Recall that the invariant density function fs [x) is defined by the equality 

fs{x) = G{S)-'e^^i^2 j^^ S{y) dyj , 
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where G {S) is the normahsing constant. 

We consider two types of tests. The first one is a direct analogue of the 
classical Cramer-von Mises and Kolmogorov-Smirnov tests based on empiri- 
cal distribution and density functions and the second follows the considered 
above (small noise) construction of tests. 

The invariant distribution function Fs {x) and this density function can 
be estimated by the empirical distribution function Ft {x) and by the local 
time type estimator fx {x) defined by the equalities 



Ft{x) 



respectively. Note that both of them are unbiased: 

EsFt (x) = Fs {x) , EsfT {x) = fs (x) , 
admit the representations 



{Xt<x} dXt 



TjT (x) = 

Ct ix) = 



Fs [Xt Ax)-Fs (Xt) Fs {x) 



Vsjx 

■ Vf 



HXt>x} 



fs (Xt) 
- Fs (Xi 



dWt + o{l), 



fs (Xt) 



dWt + o(l] 



oo 



Af{0,dFiS,xf) , 
X{0,df {S,xf) . 



and are yT asymptotically normal (as T ■ 
r^T{x) ^ y/T {Ft{x) - Fs{x)) 
Ct {x) = Vf (/t (x) - fs (x)) 
Let us fix a simple (basic) hypothesis 

J^o : S{x) = (x) . 

Then to test this hypothesis we can use these estimators for construction of 
the Cramer-von Mises and Kolmogorov-Smirnov type test statistics 



W^ = T 



D 



T — sup 

X 



Ft{x)-FsAx)\ dFsAx) 
Ft (x) - Fs^ (x) 
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and 



T 



sup 



It {x) - fs, (x) 
It (x) - fs^ (x) 



dFs^ {x) , 



respectively. Unfortunately, all these statistics are not distribution-free even 
asymptotically and the choice of the corresponding thresholds for the tests 
is much more complicated. Indeed, it was shown that the random functions 
ivrix) ,xeR) and {Ct{x) ,x & R) converge in the space ("^o,®) (of con- 
tinuous functions decreasing to zero at infinity) to the zero mean Gaussian 
processes (ry {x) ,x E R) and (C {x) ,x E R) respectively with the covariance 
functions (we omit the index of functions fs^ {x) and Fs^ {x) below) 



[V (x) T] (y)] 

[F{iAx)-F (0 F (x)] [F (e A y) - F (Q F (y)] 



Rf {x, y) 

Rf{x,y)=Y.sA(:{x)C 
= 4/ (x) / {y) 



Here ^ is a random variable with the distribution function Fs^ {x) . Of course, 

dp {S, xf = E5 [ri [xf] , df {S, xf = E5 [C [xf] . 

Using this weak convergence it is shown that these statistics converge in 
distribution (under hypothesis) to the following limits (as T — > 00) 



sup |?7 {X) 

X 



J —00 

/oo 
C{xf dF5,(x), rV2^T^SUp|C(x)|. 
-00 ^ 

The conditions and the proofs of all these properties can be found in Kutoy- 
ants (2004), where essentially different statistical problems were studied, but 
the calculus are quite close to what we need here. 

Note that the Kolmogorov-Smirnov test for ergodic diffusion was studied 
in Fournie (1992) (see as well Fournie and Kutoyants (1993) for further de- 
tails), and the weak convergence of the process ?7t (•) was obtained in Negri 
(1998). 
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The Cramer-von Mises and Kolmogorov-Smirnov type tests based on 
these statistics are 

with appropriate constants. 

The contiguous alternatives can be introduced by the following way 



S (x) = (x) + 



h (x) 

w 



Then we obtain for the Cramer-von Mises statistics the hmits (see, Kutoyants 
(2004)) 



j 2E5. [l{e<.} - Fs^ {x^)] 2^ h {s) ds^ + r/ {x) 



2fsAx) / h{s) ds + C{x) 
00 L J i 



dFs, (x) . 



Note that the transformation Yt = Fs^ (Xt) simphfies the writing, because 
the diffusion process Yt satisfies the differential equation 

dYt = fs^ (Xt) [2S. (Xt) dt + dWt] , Yo = Fs^ (Xq) 

with reflecting bounds in and 1 and (under hypothesis) has uniform on 
[0, 1] invariant distribution. Therefore, 



/ V{sfds, T^/^Dt^ sup \V{s)\, 

Jo 0<s<l 



but the covariance structure of the Gaussian process {V (s) , < s < 1} can 
be quite complicated. 

To obtain asymptotically distribution-free Cramer-von Mises type test we 
can use another statistic, which is similar to that of the preceding section. 
Let us introduce 



-I 2 



Xt-Xo- I (Xy) dv 
10 



dt. 



Then we have immediately (under hypothesis) 



Wfdt 



fw{sf 
Jo 



ds, 



where we put t = sT and W (s) — T ^^'^Wst- Under alternative we have 



Wn 



TJo 



dv 



dt 



-\ 2 



.VT TtJo 
The stochastic process Xt is ergodic, hence 
1 



dt. 



^ I h{X,)dv^^sMO^ r h(x)fsAx)dx = Ph 

^ Jo J-oo 

as t ^ oo. It can be shown (see section 2.3 in Kutoyants (2004), where we 
have the similar calculus in another problem) that 



[Phs + W{s)fds. 



Therefore the power function of the test ip {^^) = ^^w^>cc} converges 
to the function 



f3^{ph)^pQ^ [Phs + W{s)]''ds 



Using standard calcuhis we can show that for the corresponding Kolmo- 
gorov-Smirnov type test the limit will be 

(3^ (ph) = P f sup \phS + W(s)\> Cc,] . 

Vo<s<l / 

These two hmit power ftmctions are the same as in the next section de- 
voted to self-exciting alternatives of the Poisson process. We calculate these 
functions with the help of simulations in Section 5 below. 

Note that if the diffusion process is 



dXt = S {Xt) dt + a {Xt) dWt, Xo, 0<t<T, 
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but the functions S {■) and cr (•) are such that the process is ergodic then we 
introduce the statistics 



1 



T r 



Xt-Xo- d^; 



T^Es^aiCY] Jo 
Here ^ is random variable with the invariant density function 



dt. 



fs, (x) 



exp < 2 



This statistic under hypothesis is equal to 



S* jy) 
(^{yf 



dy}. 



1 



T 



T^Es. [a{0'] Jo 
1 



<j{X,)dW, 



dt 



TEs^ [a (01 Jo 



VfJo 



a (X,) dW, 



dt. 



The stochastic integral by the central limit theorem is asymptotically normal 



I ^ / a{X,)dW,^N{Q,l) 

^tEs^[a{if] Jo 



and moreover it can be shown that the vector of such integrals converges in 
distribution to the Wiener process 

(vsiT, Vs.t) =^ {W (si) ,...,W (Sk)) 

for any finite collection of < Si < S2 < • • • < -Sfc ^ 1- Therefore, under mild 
regularity conditions it can be proved that 



fw{sf 
Jo 



ds. 



The power function has the same limit, 



P^{Ph)^p(^j\phS + W{s)fds 



> Cr 
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but with 

Ph = / =■ 

The similar consideration can be done for the Kolmogorov-Smirnov type test 
too. 

We see that both tests can not distinguish the alternatives with h (•) 
such that 'Es^h{^) — 0. Note that for ergodic processes usually we have 
BsS (0 = and E^^^^/^ [5* (^) + T-'^/^h (0] = with corresponding ran- 
dom variables ^, but this does not imply 'Es^h (^) = 0. 



4 Poisson and self-exciting processes 

Poisson process is one of the simplest point processes and before taking any 
other model it is useful first of all to check the hypothesis the observed se- 
quence of events, say, < ti, . . . ,tN < T corresponds to a Poisson process. It 
is natural in many problems to suppose that this Poisson process is periodic 
of known period. For example, many daily events, signal transmission in op- 
tical communication, season variations etc. Another model of point processes 
as well frequently used is self-exciting stationary point process introduced in 
Hawkes (1972). As any stationary process it can as well describe the periodic 
changes due to the particular form of its spectral density. 

Recall that for the Poisson process X^, t > of intensity function S (t) ,t > 
we have {Xt is the counting process) 

P {Xt -Xs = k} = {k\)-' (A (t) - A (s))' exp {A (s) - A (t)} , 

where we suppose that s < t and put 

A(i) = [ S{v) dv. 
Jo 

The self-exciting process Xt,t>0 admits the representation 

Xt^ [ Sis,X) ds + TTt, 
Jo 

where ttj, t > is local martingale and the intensity function 

S{t,X)^S+ f g{t-s) dXs^S+Y.3(J-ti)- 
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It is supposed that 



POO 

/ 9 it) 

Jo 



dt < 1. 



Under this condition the self-exciting process is a stationary point process 
with the rate 

S 



1-p 



and the spectral density 



/(A) 



27r|l -^(A)! 



2 ' 



POO 

G (A) = / e'^'g (t) dt 
Jo 



(see Hawkes (1972) or Daley and Vere- Jones (2003) for details). 

We consider two problems: Poisson against another Poisson and Poisson 
against a close self-exciting point process. The first one is to test the simple 
(basic) hypothesis 



s (t) = s4t), t>o 



where S^, (t) is known periodic function of period r, against the composite 
alternative 

Mi : S{t)y^S4t), t>0, 

but S (t) is always r-periodic. 

Let us denote Xj (t) — — j — 1, . . . ,n, suppose that 

T — nr and put 

K{t)^l-J2Xj{t). 

1 1 



The corresponding goodness-of-fit tests of Cramer-von Mises and Kolmo- 
gorov-Smirnov type can be based on the statistics 



A, (r) ' n / An {t) - A, (t) dA, (t) , 



-1/2 



Dn = A* (r) -'- sup 

0<t<T 



A„ (t) - A, (t) 



It can be shown that 



WisYds, 



n Dr, 



sup |H^(s) 

0<s<l 
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where {W (s) , < s < 1} is a Wiener process (see Kutoyants (1998)). Hence 
these statistics are asymptotically distribution-free and the tests 



L{vy2>cc}, 



with the constants c^, da taken from the equations (jl]), are of asymptotic size 



a. 



Let us describe the close contiguous alternatives which reduce asymptot- 
ically this problem to signal in white Gaussian noise model (jSj). We put 



A (t) = A, (t) + 



a/^A^, (r) Jo 



A* r) 



Here h (■) is an arbitrary function defining the alternative. Then if A [t] 
satisfies this equality we have the convergence 



1 r rs 



h {v) dv + W (s) 



ds. 



This convergence describes the power function of the Cramer-von Mises 
type test under these alternatives. 

The second problem is to test the hypothesis 

: S{t) = S,, t>Q 

against nonparametric close (contiguous) alternative 



1 

S{t) = S, + ^ h{t-s)dXu t>0, 
vT Jo 



We consider the alternatives with the functions h{-) > having compact 
support and bounded. 

We have A.,, (t) = S^, t and for some fixed r > we can construct the 
same statistics 



n 



An (t) - S^t dt, Dn = (S** r) sup 

-I 0<t<T 



An (t) - t 



Of course, they have the same limits under hypothesis 



[ W{sfds, y/EDn^ sup \W{s)\. 

Jo 0<s<l 
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To describe their behaviour under any fixed alternative h (•) we have to 
find the hmit distribution of the vector 

1 " 

Wn = (Wn (tl) ,...,Wn {tk)) , Wn {tl) = [Xj {ti) - SJi] , 

VS*T n 

where < t/ < r. We know that this vector under hypothesis is asymptoti- 
cally normal 

>Co{w„}^Ar(0,R) 

with covariance matrix 

R- = (-R/m)jkxjt ' Rim — 'T ^ min (t;, i^) . 

Moreover, it was shown in Dachian and Kutoyants (2006) that for such al- 
ternatives the likelihood ratio is locally asymptotically normal, i.e., the like- 
lihood ratio admits the representation 

Z„ (h) = exp I A„ {h, X^)-h (h) + r„ {h, X^) | 

where 

A„ {h, X^) = / / h{t-s) dX, [dXt - S^dt] , 



l{h) 

and 



roo / roc 

J h{tfdt + sjj h{t)dt 



A„ (h, X") ^X{0,1 (h)) , r„ (h, X") ^ 0. (7) 

To use the Third Le Cam's Lemma we describe the hmit behaviour of the 
vector (A„ {h, X") , w„). For the covariance Q = {Qim) ,l,m — 0,1, ... ,k oi 
this vector we have 

EoA„ {h, X") = 0, Qoo = EoA„ {h, X")' = I (/i) (1 + o (1)) . 

Further, let us denote dvrt = dXt — S^dt and H (t) = h{t — s) dXg, then 
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we can write 



Qo^ = Eo [A„(/i,X")ti;„ (ti)] 



, = 1 -'rij-l) 
r(j-l)+ti 



i=l 



-{i-l)+ti 



T{i-l) 



, EoH {t) dt ^ ^ ^/K h{t)dt {l + o{l)), 

nrVS* ^ Jr{j-i) T 



EqH {t) = S, h{t~ s) ds = S, h (s) ds 



because 



for the large values of t (such that [0, t] covers the support of h{-)). 
Therefore, if we denote 



h 



/•oo 

- / his) 
Jo 



ds 



then 



Qoi — Qio — — \/sl h. 

T 



The proof of the Theorem 1 in Dachian and Kutoyants (2006) can be 
applied to the linear combination of A„ (h, X") and Wn (ti) , . . . ,Wn (tk) and 
this yields the asymptotic normality 



£o(An(/i,^"),w,) ^AA(0,Q). 



Hence by the Third Lemma of Le Cam we obtain the asymptotic normality 
of the vector 



[v/nj =^ \ W (si) + si V 5'* /i, . . . , Vl/^ (sfc) + Sfc V 5"* /ij , 

where we put ti — tsi. This weak convergence together with the estimates 
hke 

E/j \Wn (ti) - Wn it2)f < C \ti - tgl 

provides the convergence (under alternative) 



S,hs + W{s) 
17 



ds. 



We see that the hmit experiment is of the type 

dYs = hds + dW{s), Yo = 0, < s < 1. 

The power jS^t/jn, h) of the Cramer- von Mises type test ipn^X^) = 1{vk2>cq} 
is a function of the real parameter ph = \/S^ h 

(3 {Wn. h) = p(^j^ [Phs + W {s)f ds >c^^+o (1) = {ph) + (1) . 

Using the arguments of Lemma 6.2 in Kutoyants (1998) it can be shown 
that for the Kolmogorov-Smirnov type test we have the convergence 

A/nD„ =^ sup \phS + W (s) I . 

0<s<l 

The hmit power function is 

134, (ph) = P ( sup \phS + W {s)\ > da) . 

Vo<s<l / 

These two hmit power functions will be obtained by simulation in the 
next section. 



5 Simulation 

First, we present the simulation of the thresholds Ca and da of our Cramer- 
von Mises and Kolmogorov-Smirnov type tests. Since these thresholds are 
given by the equations (j4]), we obtain them by simulating 10^ trajectories of 
a Wiener process on [0,1] and calculating empirical 1 — a quantiles of the 
statistics 



W'= / Wisyds and D= sup \W {s)\ 

Jo 0<s<l 

respectively. Note that the distribution of W'^ coincides with the distribution 
of the quadratic form 

oo .2 

^' = E7T^' i.i.d. ~ AT (0,1) 
fc=i y^'^) 
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and both distributions arc extensively studied (see (1.9.4(1)) and (1.15.4) in 
Borodin and Salmienen (2002)). The analytical expressions are quite com- 
plicated and we would like to compare by simulation Cq, and da with the real 
(finite time) thresholds giving the tests of exact size a, that is and 
given by equations 

P {W^ >cl} = a and P {v^D„ > dl} = a 

respectively. We choose S* = 1 and obtain and d^ by simulating 10'' tra- 
jectories of a Poisson process of intensity 1 on [0,T] and calculating empirical 
1 — a quantiles of the statistics and y/nDn- The thresholds simulated for 
T = 10, T = 100 and for the limiting case are presented in Fig. 1. The lower 
curves correspond to the Cramer-von Mises type test, and the upper ones 
to the Kolmogorov-Smirnov type test. As we can see, for T = 100 the real 
thresholds are already indistinguishable from the limiting ones, especially in 
the case of the Cramer-von Mises type test. 




0,2 0,4 0,6 0,8 1 

alpha 

Fig. 1: Threshold choice 

It is interesting to compare the asymptotics of the Cramer-von Mises and 
Kolmogorov-Smirnov type tests with the locally asymptotically uniformly 
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most powerful (LAUMP) test 



nr 



proposed for this problem in Dachian and Kutoyants (2006). Here Za is 1 — a 
quantile of the standard Gaussian law, P > Za) — a, C ~ ■A/'(0, 1). The 
limit power function of is 

In Fig. 2 we compare the limit power functions [3^ (p) , jS^f, (p) and (3^{p)- 
The last one can clearly be calculated directly, and the first two are obtained 
by simulating 10^ trajectories of a Wiener process on [0,1] and calculating 
empirical frequencies of the events 

< / [ps + W {s)f ds > Ca\ and < sup \ps + W {s)\ > da 
{Jo J lo<s<i 

respectively. 




Fig. 2: Limit power functions 
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The simulation shows the exact (quantitative) comparison of the Hmit 
power functions. We see that the power of LAUMP test is higher that the 
two others and this is of course evident. We see also that the Kolmogorov- 
Smirnov type test is more powerful that the Cramer- von Mises type test. 
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